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ABSTRACT 

The role of turbulence in a spherically symmetric accreting system has been 
studied on very large spatial scales of the system. This is also a highly subsonic 
flow region and here the accreting fluid has been treated as nearly incompressible. 
It has been shown here that the coupling of the mean flow and the turbulent 
fluctuations, gives rise to a scaling relation for an effective "turbulent viscosity" . 
This in turn leads to a dynamic scaling for sound propagation in the accretion 
process. As a consequence of this scaling, the sonic horizon of the transonic 
inflow solution is shifted inwards, in comparison with the inviscid flow. 

Subject headings: accretion, accretion disks — hydrodynamics — methods: ana- 
lytical — turbulence 



1. Introduction 

The purpose of this work has been to study the dynamic scaling behaviour of the 
coefficients of viscosity arising out of turbulence in a spherically symmetric accreting system, 
and how such scaling behaviour leads to a scale dependence for the speed of sound as well. In 
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its turn this will be shown to have an important bearing on the sonic point of the transonic 
inflow solution, since it is with the speed of sound that both the bulk velocity of the flow 
and the sonic point are scaled. A study of this kind should be useful in addressing recent 
observational discrepancies, for which the classical (and inviscid) Bondi theory has proved 
somewhat inadequate. 

In accretion studies, turbulence is of great relevance, since in almost all cases of physi- 
cal interest, the accreting astrophysical fluid is in a turbulent state (Choudhuri 1999). It is 
presently a well established fact that in the two dimensional case of a thin accretion disc, vis- 
cous shearing between two differentially rotating adjacent layers, accomplishes the outward 
transport of angular momentum and effectively facilitates the infall of matter (Pringle 1981; 
Frank et al. 1992). However, in this situation, ordinary molecular viscosity has been known 
to be quite an inadequate mechanism to explain the rate of the transport process. On the 
other hand, it is to a very high value of the Reynold's number that the flow corresponds, 
and as such the flow is widely acknowledged to be turbulent (Frank et al. 1992; Balbus & 
Hawley 1998). In such a situation, turbulence — as quantitatively characterized by a "tur- 
bulent viscosity" in the Navier-Stokes' equation — becomes a prime candidate for a physical 
mechanism that brings about an enhanced outward transport of angular momentum. The 
very well known a prescription of Shakura & Sunyaev (1973) is based on this principle. 

As opposed to the facilitating role that it plays in a rotationally accreting flow, viscosity 
- even presumably "turbulent viscosity" — affects the more paradigmatic spherically sym- 
metric accreting flow, somewhat differently. In the latter case, it has been seen that the role 
of viscosity is actually directed towards inhibiting the process of gravity driven infall, and 
in doing so, viscosity also sets up a limiting length scale for the effectiveness of gravity - 
the "viscous shielding radius" (Ray 2003). However, as in the thin disc system, even in the 
spherically symmetric case, molecular viscosity would be far too weak a mechanism to bring 
about a significant quantitative impact. It would then be well worth investigating into the 
question of how significantly would spherically symmetric accretion be affected by a large 
and scale dependent turbulent viscosity. Following the qualitative insights obtained with the 
introduction of molecular viscosity in the governing hydrodynamical equations, it is possible 
as a matter of standard practice to study both the qualitative and the quantitative extent of 
the influence of an effective turbulent viscosity on the hydrodynamical processes. The main 
purpose here would be to show that the turbulent fluctuations of the interstellar medium 
are capable of renormalizing on large length scales, the small molecular viscosity given in 
the Navier-Stokes equation. This renormalized effective viscosity, as pictured by Heisenberg 
(1948) in his theory of turbulence, can very well be instrumental in setting a noticeable 
limiting length scale on the effectiveness of gravity to drive the accretion process. Indeed, 
the renormalizing of the viscosity would be robust enough to make viscosity be comparable 
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with pressure, on the same scale of length. 

Related to this contention, one example that may be cited is that of the resistive role 
against gravity that turbulence plays in another spherically symmetric system of astrophys- 
ical interest — that of the self-gravity driven Jeans collapse of a gas cloud without angular 
momentum. Studies carried out by Bonazzola et al. (1987, 1992) have shown that turbulence 
acts as a stabilizing agent against a self-gravity driven collapse. A renormalization approach 
has shown that a renormalized turbulent pressure acts against gravity. In addition to this, 
what is being contended in this work is that a scaled-up renormalized viscosity also enfeebles 
the influence of gravity. Both the pressure term and the viscous term derive from the stress 
tensor in the Navier-Stokes equation (Tennekes & Lumley 1972). For the spherically sym- 
metric case, the contribution comes only from one diagonal (hence, isotropic) element of the 
stress tensor. In such a situation, the renormalized pressure and the viscosity terms would 
both manifest themselves through the same physical effect, and the physical contents of the 
arguments presented by Bonazzola et al. (1987, 1992) for the self-gravity driven spherical 
collapse of a gas cloud, would match those in this study of turbulent spherically symmetric 
accretion. Having noted this point, it would also be instructive to have an understanding 
of the difference between the two physical cases being compared here. Whereas Bonazzola 
et al. (1992) have studied the response to large scale density perturbations on a stationary 
turbulent solution in a self-gravity driven system, what is being studied in this work, is the 
influence of spontaneous fluctuations on the mean stationary solution of a system, in which 
gravity comes into play through an external accretor. An analysis of the latter nature is all 
the more contextual with regard to accretion studies, because spherically symmetric accre- 
tion is exemplified by the infall of interstellar matter on to an isolated accretor, and it has 
been well recognized that the interstellar medium displays turbulent behaviour (Jokipii & 
Lerche 1969; Jokipii et al. 1969; Lee & Jokipii 1976). 

While dwelling on this matter, it would be important to mention that certain previous 
studies in spherically symmetric accretion on to a black hole, have in fact quantitatively 
accounted for the physical role of turbulence in very efficiently converting gravitational energy 
to radiation. In the works of Meszaros (1975) and Meszaros & Silk (1977), it has been argued 
that for spherical accretion on to a massive black hole, turbulent dissipation would be one 
of the factors which would result in the luminosity of the system being enhanced by quite a 
few orders of magnitude — indeed to such an extent that the spherically symmetric system 
could be compared with disc models as an X-ray source. 
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2. The equations of turbulent spherical accretion 

The effect of turbulent fluctuations has been studied here on very large spatial (and 
therefore highly subsonic scales) of the spherically symmetric accreting system. The physical 
effects of turbulence are appreciably manifested on these scales. Since all physically feasible 
flow solutions have to pass through the subsonic flow region, the turbulent fluctuations here 
must have a significant influence on the flow. And more to the point, on these subsonic 
scales, the flow could be studied in the nearly incompressible regime. 

Turbulence is an attribute of the fluid flow (Tennekes & Lumley 1972), while molecular 
viscosity is an intrinsic physical property of the fluid. And yet the two can be very closely 
related to each other through the Navier-Stokes equation (Frisch 1999), which, as one of the 
governing equations of the flow, is given by 

+ (v-V)v + + r = z/V 2 v + /iV(V-v) (1) 



dt p r 2 

where v and p are the two kinematic coefficients of viscosity. The pressure P is related 
to the density through a general polytropic equation of state P = kp 1 . Here 7 is the 
polytropic exponent with an admissible range given byl<7<5/3 — these restrictions 
having been imposed by the isothermal limit and the adiabatic limit respectively. The flow 
is also governed by the continuity equation, which is given by 

f + V-( P v) = (2) 

The total velocity and density fields are written as v = v + u and p = Po + $P, i n which 
Vo and po, which are functions of the radial coordinate only, are the mean velocity and 
density profiles for the spherically symmetric transonic flow, while u and Sp generally are 
time-dependent and three-dimensional random fluctuations about the transonic solution. 
The implicit understanding here is that in an accreting system naturally evolving in real 
time, the transonic solution is accorded primacy over all possible other stationary solutions 
in both the inviscid (Bondi 1952; Garlick 1979; Ray & Bhattacharjee 2002) and viscous 
regimes (Axford & Newman 1967). Under the assumption that cross-correlations of the 
density and the velocity fluctuations would be negligible, i.e. (V-(u5p)) = 0, the average 
(and steady) solutions would be obtained as 

V-(poVo) = (3) 

and 

(vo-V)vo + ((u-V)u) + ^j^Po + — r = z/V 2 v + pV(V-v ) (4) 
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In the nearly incompressible regime, only the first order term in the expansion of the 
density fluctuations about the mean density need be retained. The fluctuating density and 
velocity fields are therefore seen to satisfy, 



d 

—Sp + V-(pou) + V-(v 5p) + V-(u5p) 







(5) 



and 



du . _ 
^ + (u.V)u + V 



,5p 



z/V 2 u + pV(V-u) - (u-V)vq + (vo-V)u - ((u-V)u) (6) 



respectively, with c s being the steady value of the speed of sound, which is related to the 
mean density by cj. = •ykp^ 1 . 

In the very much subsonic region of the flow, the variation of the mean density po ma Y 
be neglected, since in this region the mean density very closely assumes an ambient value, 
which is a constant. Under this approximation, equation (3) gives the relation V-v = 0. 
Furthermore, on these scales the continuity equation also governs the asymptotic behaviour 
of the mean velocity, which implies that the variation of the mean velocity, at its most rapid, 
is given by v ~ r~ 2 (Petterson et al. 1980; Chakrabarti 1990). On the other hand (under 
these asymptotic conditions) the turbulent velocity fluctuations are much greater than the 
mean velocity itself, and in fact are of the order of the speed of sound. Hence on these large 
length scales, ignoring all terms involving the mean velocity and the gradient of the mean 
density and its fluctuations, it would be meaningful to retain only the primary signature of 
a compressible flow, namely V-u ^ 0, whence equation (5) simplifies to 



1 d 

— ^p + V-u = 

Po ot 



(7) 



which is an expression that has found quite regular mention in the study of a nearly incom- 
pressible fluid flow with random fluctuations (Staroselsky et al. 1990; Bhattacharjee 1993). 

At this stage it should be important to be assured of the consistency in neglecting 
the higher powers of Sp/po. Under the chosen working approximations, the terms in the 
left hand side of equation (6) can be written as u a , updpu a and c 2 s d a (5p/ po) respectively. 
If u a and c s are to scale as L € , then the time t scales as L 1 ^ 6 , while (Sp/po), of course, 
remains independent of any scaling. Here e is arbitrary, but anticipating that a one-loop 
calculation will yield a positive value for e, all nonlinearities involving Sp/po (with p being 
asymptotically a constant) may be ignored in favour of Updpu a . This, arguably, should 
suffice for a study of the scaling dependence in the flow. If the resulting calculations lead to 
a positive e, the adopted procedure would be justified and would be consistent with itself. 
That this is precisely what happens, will be demonstrated in the following sections. 



- 6- 



Equation (6) is likewise simplified, and closed with the help of equation (7), to give 

— + (u-V)u = C 2v(-J (V-u) + z/V 2 u + ^V(V-u) + f (8) 

in which f = — [(u-V)vo + (vo-V)u — ((u-V)u)]. The primary complication at this stage is 
in this term f , which couples the fluctuating flow to the mean flow. Gravity of the central 
accretor maintains the mean flow, from which energy is transferred to the fluctuating flow, 
through its coupling with the mean flow. So effectively what happens is that the turbulent 
fluctuations are sustained by gravitation, via the nonlinear coupling in the term f . Various 
approximations in the theory of turbulence have involved a modelling of this nature of energy 
input to the turbulent flow. Prandtl's mixing length theory is one of the most well known 
(Faber 1995; Choudhuri 1999). A more recent point of view treats this force as an as yet 
unspecified force external to the turbulent flow (Forster et al. 1977; De Dominicis & Martin 
1979). Its dependence on the random field u, makes it random and hence the modelling 
endows f with random properties. Even for this accretion problem it would therefore be 
quite possible to conceive of a randomly forced turbulent flow described by (for the nearly 
incompressible flow that is being studied here) 

d t Ui + (ujdj)ui = c^diid^djUj) + udjdjUi + fid^djUj) + (9) 

in which, for the Gaussian forcing, the correlation function is specified as 

(f t (r,t)f 3 (r,t)) = 5 l3 C (\r-r>\)5(t-t>) (10) 

These two equations (9) and (10) will be necessary to develop a dynamic scaling theory for 
the turbulent spherically symmetric flow. 



3. Dynamic scaling for turbulent spherical accretion 

To carry out a dynamic scaling analysis with the help of equations (9) and (10), it would 
be convenient to work in Fourier transform space. This would necessitate writing 



u 



(r,t) = j^Ju t (k,uj)e^-^d 3 kduj (11) 



in terms of which equation (9) becomes 

1 kk- v-^ 

—vjj + uk 2 )5ij + fikikj Uj — (? s -r-^-Uj — f% — i y^PjUj(k — p, u/)iij(p, uo — uo') (12) 

p,<y 
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The technique that has been adopted here is to expand the velocity field as Ui = uf^ + 
u*p + u\ 2 ^ + . . . , in which uf^ is the solution in the absence of the nonlinear term. The 
subsequent terms are the effect of the nonlinear term in equation (12). The lowest order 
solution can then be written as 

<4 0) = GV> fj (13) 

in which 



,(0) ] 1 = ( - ico + z/A; 2 ) 5ij - - /i) hkj (14) 
The first-order correction, uf\ satisfies 

2 

( - icu + vk 2 ^ - - vtykikj] uf ] = -iJ2Pj u f\ k ~ P> ^')4 0) (P, w - c/) (15) 



and its solution is given by 

uf ] = -iG\f(k,u) ^p fc 4 0) (k - p,u')uf\p,u - J) (16) 



As has been stressed by Heisenberg (1948), the momentum transfer term, given in the 
right hand side of equation (12), gives rise to an effective turbulent shear viscosity (the 
eddy viscosity) for an incompressible flow. This is the physical content of all subsequent 
theories — the different kinds of renormalized perturbation expansion (McComb 1990), the 
renormalization group (McComb 1990), the self-consistent mode coupling (McComb 1990) 
and the very recent Lagrangian picture approach (L'vov & Procaccia 1995a,b). In this 
compressible case, it will be easy to see that the right hand side of equation (12) will give 
rise to the effective shear viscosity, bulk viscosity and speed of sound. The simplest way of 
arriving at this result is to examine the average value of the right hand side of equation (12), 
averaged over the distribution of the random force /j. This is done perturbatively. In what 
follows in this section, the salient results of the perturbative analysis have been brought 
forth. The details of the calculations have been presented in the Appendix. 

The averaging of the nonlinear term in equation (12) will lead to its equivalent linearized 
representation (see the Appendix), given by 

(-i^2pjUj{k-p,u')ui(p,u-u')) = -a^\k,u)uf\k,uj) (17) 

p,w' 

where 

4 0) (k, u) = 2j2Pjh [Gf(k - p, cv')C^(p, u - J) + 4 0) (p, u - o/)<7g>(k - p, a/)] (18) 

p,uj' 
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and in which 



(19) 



Clearly, is the zeroth-order correlation function for the velocity field. The nonlinear 
term in the equation of motion now has the structure <r^(k, o;)u|°^(k, uj) in this lowest 
order of the perturbation theory. The coefficient can clearly be identified as making a 
contribution to the two coefficients of viscosity and the speed of sound, by comparing with 
the linear term in the equation of motion shown in equation (12). The that has been 
obtained is called the self energy and constitutes a dressing of the bare coefficients. This 
is exactly in conformity with all the different ways of doing perturbation theory (McComb 
1990). The step beyond perturbation theory goes to say that as the higher order terms are 
considered, will be converted to the full self energy an. It must be emphasized here that 
conversion to the full self energy will not differently affect the scaling arguments that will 
be developed here on the basis of the lowest order in the perturbation theory. 

To make any further progress, it would be instructive to examine the structure of 
cj^(k, a;)u[ ^(k, uj). By comparing with the form of [G^]" 1 in equation (14), it is possi- 
ble to write 



where <j\ and a 2 are the frequency and momentum dependent components of the self 
energy tensor which must have the structure shown in equation (20) from the isotropy of 
space. Evidently, cr^(k, uj = 0) dresses the shear viscosity, while cr^k, uj) dresses the bulk 
viscosity and the speed of sound. To have any information about the dressing of the speed 
of sound, the (iuj)~ l part would have to be extracted from a^(k, u) and the rest of the 
integral would have to be evaluated at uj = 0, to yield the dressed bulk viscosity. 

The renormalization of z/, fi and c s converts them into the renormalized quantitites u, ji 
and c s respectively. To have any idea of how the two coefficients of viscosity and the speed of 
sound get renormalized, the Green's function would have to be written out by inversion of the 
matrix implied by [G^] _1 in equation (14). Substitution of the unrenormalized quantities 
in the Green's function by the renormalized ones (see the Appendix), will then give the fully 
dressed Green's function as 



alfiKu;) = 2j2PmK G^](k-p,uj')C^(p,uj-uj') 




(20) 



1 




(21) 



iuj + uk 2 + k 2 (fi — c 2 /ioj) 



9 



The poles of the Green's function, occurring at u — —ivk 2 and the roots of uo 2 = 
c 2 k 2 — iuk 2 \v + jl) , would deliver a dispersion relation for uj. It is satisfying to note here that 
the second relation (the quadratic in u) is identical to the one obtained by Bonazzola et al. 
(1992), barring a term arising from the self-gravity of the system that they were studying. 
In the long wavelength limit (k small), which contains the interesting features about the 
scaling behaviour, the quadratic in u can be approximated as 



The renormalized quantities will have a power law dependence on k. If dynamic scaling 
is to be invoked, then the frequency will be proportional to some definite power of k, which 
means that all the terms in the right hand side of equation (22) must scale in the same 
way. In physical terms this would mean that both the propagating term and the dissipative 
term in equation (22), would be comparably effective on the same scale. If a power law 
were to be written in the form v oc k~ y , then it will clearly also indicate that ji oc k~ y 
and c s oc k 1 ~ y '. The main concern will be to set a value for y. To make any progress in 
that direction, the forcing term would have to be specified, which actually would imply 
specifying the correlation function C (|r — r'|) in equation (10). Since this term dominates 
at large distance, a scaling form C (|r — r'|) oc |r — r'\ a may be assumed. This would then 
transform in the momentum space as Co(k) ~ k~ ( - D+a \ where D is the dimensionality of the 
space. At this point a significant departure from Bonazzola et al. (1992) is being made, by 
suggesting that both the renormalized af^ and as given by equation (20), would have 
to be treated as comparable with each other. Consequently, information on the scaling of 
the sound velocity (arising from the pressure term) can be had from the scaling of the shear 
viscosity, which is clearly a dissipative effect. 

The fully dressed self-consistent form of equation (20) can now be written down as an 
integral (see the Appendix), given by 



where C ni = G nr C G* ri . 

For an incompressible flow, the Kolmogorov spectrum requires that a — 0, to charac- 
terize the nature of the transfer of energy between the mean flow and the fluctuating flow. 
Since this transfer characteristic should be independent of the speed of sound, it would be 
possible to write a = for the case of near incompressibility being discussed here. The left 
hand side of equation (23) then scales as k 2 ~ y , while the right hand side scales as k 2y ~ 2 . For 
the scaling properties of both sides of equation (23) to agree, it should be necessary to set 



uj = ±c s k k 2 (v + Jx) 



(22) 




(23) 
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2 — y = 2y — 2, which will yield y — 4/3. This will then establish the result 

v~k^l\ ii~k~ A /\ ~c s ~k-V 3 (24) 

which, it may be mentioned at this point, is identical to the scaling relation obtained by 
Staroselsky et al. (1990) for the case of a randomly stirred compressible fluid. It must also 
be emphasized here that simple dimensional arguments would not entirely suffice. Indeed, 
Staroselsky et al. (1990) make this point amply clear by saying that the renormalization of 
the speed of sound is essential to understanding the physics of compressible flows, since the 
appearance of the speed of sound as a dimensional parameter, makes simple dimensional 
considerations invalid. 

Of immediate interest would be the scaling behaviour of the speed of sound, which 
in terms of the radial distance may be written as c s ~ r 1 / 3 . The steady state solution of 
the continuity equation (with p written in terms of c s ), gives a dependence for the steady 
velocity of the flow v , which goes as (Chakrabarti 1990) 

v ~ r- 2 c~ 2n (25) 

where n = (7 — l)" 1 is the polytropic index, whose admissible range of values for inflow 
solutions is given by3/2<n<oo (Chakrabarti 1990). Using the renormalized speed of 
sound and its associated scaling relation in equation (25), will give a scaling behaviour for 
the steady flow velocity as 

Vo „ r -2(l+n/3) (26) 

from which it is quite evident that regardless of the value of n, on large length scales, the 
steady flow velocity would die out — a fact that is in conformity with the boundary condition 
of the flow. The result given by equation (26) highlights another very interesting issue. It has 
been discussed earlier that on large length scales, the mean flow is limited by the equation 
of continuity, and therefore its variation is given by Vq ~ r~ 2 . This is a result that is easily 
derived from the classical and inviscid Bondi flow (Petterson et al. 1980; Chakrabarti 1990). 
What equation (26) indicates is that turbulent fluctuations, sustaining themselves at the 
expense of the mean flow, detracts even further from the r~ 2 scaling law for the mean flow 
velocity — something that, from considerations of energy dissipation, can be qualitatively 
intuited about the influence of turbulence on the mean flow. 

The results in equations (24) and (26) also lead to the conclusion that in the renormalized 
situation, there would be a scale dependence for the position of the sonic horizon as well. 
For large length scales, i.e. concomitantly for a large effective turbulent viscosity, the sonic 
horizon would be shifted inwards. This happens because, seen on a large length scale, an 
enhanced scale dependent speed of sound, could only be matched by the steady flow velocity 
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deeper within the gravitational potential well. Since the flow has to pass through the subsonic 
region in any case, this effect of subsonic turbulence in shifting the sonic point inwards, is 
also seen to have a bearing on the transonicity of the inflow solution. This observation would 
be entirely compatible with the role of a weak molecular viscosity in inwardly shifting the 
position of the critical point of the inflow solution (Ray 2003). 

4. Concluding remarks 

It has been seen so far, how on the large length scales of a spherically symmetric accreting 
system, turbulence is capable of setting a scaling behaviour for both viscosity and the speed 
of sound. However, it need not be supposed that given the scaling relation c s ~ r 1 / 3 ) 
there would be an arbitrarily large scaling for the speed of sound on large length scales. 
This is because in spherical symmetry, turbulence itself will also play a role in limiting the 
accretion process. The physical quantity m/up (with m being the accretion flow rate) has 
the dimension of length, and this has been understood to be a viscous shielding radius, r visc 
(Ray 2003). If the value of v is enhanced by the introduction of a large and scale dependent 
kinematic viscosity, then r visc will define a noticeable spatial limit for the accretion process. 

The r 1//3 scaling behaviour for sound propagation is also apparently surprising, with its 
physical implication being that the flow is heated up more at larger radial distances. On the 
other hand, the classical Bondi theory shows that the speed of sound increases as the flow 
moves inward, i.e. the flow gets heated up more at smaller radii. The point to remember 
is that this property of classical spherical accretion is not violated by the mean flow, and 
the dressing of sound propagation that the turbulent fluctuations bring about, is manifest 
over and above the standard features that the mean flow is expected to show. The extent 
of energy dissipation that turbulence brings about is not accounted for by the Bondi theory. 
This energy dissipation shows itself as an enhanced scaling for the speed of sound (larger 
scales are more energetic in this sense), and had temperature been chosen as a dynamical 
variable, this would have shown no contradiction. A cautionary reminder that is to be 
sounded here is that all the scaling relations have been derived under the assumption of near 
incompressibility on large length scales, which is a condition that cannot be applied too far 
into the inner region of the flow, and in consequence, the r 1 / 3 scaling for sound propagation 
is not to be extended too much to small length scales either. 

It would also be instructive here to have an understanding of the dynamic scaling of 
both the speed of sound and the steady flow velocity, which could be derived on using the 
prescription for an effective viscosity forwarded by Meszaros & Silk (1977). They proposed 
a scaling behaviour for the kinematic coefficient of turbulent viscosity u t , which could be 
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conceived of as a product of a characteristic length scale l t of the turbulent cells, and the 
magnitude of their associated turbulent velocity fluctuations v t . It was assumed by Meszaros 
& Silk (1977) that for all r, the length / t would be some fraction of the radial distance r, 
while v t would be a fraction of the free fall velocity vr , which varies as r -1 / 2 . For large length 
scales, in which the density of the accreting fluid approaches its constant ambient value, this 
prescription would lead to a scaling behaviour given by v t ~ Vtk ~ r 1 / 2 . 

In this case, it would then be easy to see from the dispersion relation given by equation 
(22), that the speed of sound would be scaled by the relation c s ~ r -1 / 2 , while scaling for 
the steady flow velocity, from equation (25), would be given by vq ~ r n ~ 2 , The difficulty 
arises for n > 2, since for large length scales, v would actually increase, contrary to a 
common understanding of the boundary condition that v should decrease over large radial 
distances. This discrepancy arises because of considering the characteristic eddy velocity to 
be a fraction of the free fall velocity. Even though this looks well founded on dimensional 
principles alone, this scaling behaviour breaks down on large length scales, because on these 
scales free fall conditions do not hold. Rather, this is the region of the ambient conditions, 
where the mean flow velocity, even under inviscid conditions, varies at the most as r~ 2 , and 
therefore the velocity fluctuations would have to have a different scaling behaviour. Indeed, 
by being coupled to the mean flow, the velocity fluctuations alter the scaling behaviour for 
the mean velocity as well, as equation (26) indicates. Thus it would probably be more correct 
to suggest that within the sonic radius and close to the accretor, for a highly supersonic mean 
flow, free fall conditions can have a bearing on the velocity fluctuations. However, the extent 
of the influence of turbulence on such small scales would be a somewhat contentious issue, 
and is not within the scope of this work. 

This research has made use of NASA's Astrophysics Data System. One of the authors 
(AKR) gratefully acknowledges the support provided by the Council of Scientific and Indus- 
trial Research, Government of India, for a part of the time that was needed to carry out this 
work. 



Appendix 



It has been seen that equation (12) is in the form 



{—iu + uk 2 )5ij + fikikj Uj — c^-^-Uj — fi — i ^^PjUj(k — p, c</)-Uj(p, u — uo') (Al) 



p,u/ 
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of which, the right hand side is averaged over the distribution of the random force /j. For 
the nonlinear term, the perturbative expansion of Ui can be written for its first two terms as 

(-i^2p jUj (k - p, u')ui{p, u-u')) = (-i^pjufik - p, u')uf\p, u - u')) 

+ h E ps i*? ( k - P' w 'K (0) (p. w - J ) 

p,uj' 



+ u f (k-p,c'K (ij (p^-^)]) 

(A2) 

In the above equation, the first term on the right hand, with uf^ substituted from 
equation (13), can be written as 

<-i (k - p, u/)ii<°>(p, u-u')) = -i £p,-Gg(k - P> "O^P, * - a;') 



p,w' 



X (/ m (k - p, w')/n(p, ^ - u')) 

p,u>' 

As can easily be seen, the expression above does not produce any momentum (k) or 
frequency (u) dependent term and hence is not responsible for momentum transfer. The 
second term in equation (A2), with u^p substituted from equation (16), can be written down 

as 

(-i E Pj [ u f ] ( k - P> L0> ) u i° ) (P> w - w') + u f ( k ~ P» UJ ') ut t ) (P» u ~ UJ ')]) 

P,U)' 

= - <E ^ [ G fi ( k - P' w ') E <^ 0) ( k - p - * ^>f 0) fa w ' - w >i 0) (p» w - ^ 

p,w' q,w" 

+«f (k - P , c')4 0) (p, w - w') E (p - q. ^>! 0) (q. ^ - a/ - a/')] ) 

q,w" 

-2 E Pi[^ ) (k-p,a;0?*C^(k-p-q,OG£ ) (p,a;- t y) 



x (/ m (k - p - q, w")/„(p, w - t/))uj 0) (q, - u") 



+4 0) (p, a; - «/)q k GZ(P ~ * <0<3» ( k - p, a/) 



,(o) 



«\r»(o), 



X (/m(P - q, w")/n(k - P> ^')) M S 0) (q) W-w'-W 



The factor of 2 appears in the expression above because -u^ could be expressed in two 
ways with the help of equation (13). Using the correlation function implied by equation (10), 
will now give from the result above 
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-2 J^Pjh \Gf(k - p, u,')G^(-p, -u + u/)G£>(p, u; - u/)C (p) 



t(0). 



+GS"(p, c - u/)G£>(p - k, -V)GjS(k - p, u/)C (k - p)j u\ u > (k, W ) 
-2 (k - P, tS)C£\p, u - u/) 



(o). 



.(0)/ 



+4 0) (p, w - u/)Cg J (k - p, J) j (k, W ) 

= -4 o) (k,c) M ; o) (k,c) 

where 



y(0), 



,(°)/ 



G<? (k - p, a/)C7ff (p, a; - a;') + (p, u - a/)CT (k - p, a/) 



(0), 



y(0). 



P,U!' 



and 



In the lowest order of the perturbation theory, the nonlinear term now has an equivalent 
linearized representation given by <r-^(k, u;)u^(k, u). The that has been obtained is 
called the self energy and it serves the purpose of dressing the bare coefficients in the equation 
of motion. Considering all the higher order terms, will be converted to the full self 
energy an. The self energy can be compared with equation (14) and can be seen to make a 
contribution to the two coefficients of viscosity and the speed of sound. Seen in this way it 
can be written as 



a^(k,co) = 2j2Pmk n G^(k-p,u/)(7SW-^) 



,(0), 



+^(P^-^C^(k-p,a/) 



(0), 



k 



a^ (k, uo)Sij + a^' (k, uo 



(0), 



k 2 J 



(A3) 



in which af\k,u = 0) dresses the shear viscosity, and a^'(k,uj) dresses the bulk viscosity 
and the speed of sound. 

To understand the effect of renormalization, it would be necessary first to obtain the 
Green's function by inversion of the matrix implied by equation (14). In this way the Green's 
function is given as 

(n-c 2 Jiu) 



(0), 



G\f(k,u) 



—iuj + vk 2 



( y . _ fc.fc _ 

lJ 10 —iu + vk 2 + k 2 (fi - c 2 Jiuj) 



(A4) 
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It is important to check if the incompressible limit is to be correctly obtained. The in- 
compressible limit implies c s — > oo and in that limit it is seen that G^(k, uj) = Pjj(k)[— iu+ 
uk 2 ]' 1 , where i^-(k) = 5^ — (kikj)/k 2 , is the projection operator, which is as it should be. 

The renormalization of u, \i and c s converts them into the renormalized quantities i>, ft 
and c s respectively. In the event of the two coefficients of viscosity and the speed of sound 
getting renormalized, the fully dressed Green's function is given by 



—iu + vk 2 



(fi-c 2 s /iuj 



,y . _ /,. /,. " .. 

%3 1 3 —iu + vk 2 + k 2 (ft- c 2 /iuj) 



(A5) 



The poles of the Green's function occur at uj = —ivk 2 and the roots of uj = c s k 
iu)k 2 (u + p), and on solving the quadratic in cu, the dispersion relation is given by 



2u = -ik 2 {y + ft) ± ^Ac 2 s k 2 - k\v + fx) 2 (A6) 
The long wavelength limit (k small) yields 



uj = ±c s k- l -k 2 (v + jj) (A7) 

Dynamic scaling would imply that the frequency would be proportional to some power 
of k, and this in its turn would mean that each term in the right hand side of equation (A7) 
must scale in the same way. Assuming a power law of the form v oc k~ y , will also clearly lead 
to having ft, oc k~ y and c s oc k l ~ y . To know the value of y, the function C (|r — r'|) in equation 
(10) would have to be specified first. Assuming a scaling form given by Co(|r— r'|) oc |r — r'\ a , 
yields the corresponding transformation in the momentum space as Co(k) ~ k~^ D+a \ with 
D being the dimensionality of the space. 

The fully dressed self-consistent form of equation (A3) can be set in an integral form 
given by (on noting that the two summations in the right hand side of equation (A3) are 
quite identical) 

<Tij(k,uj) = 4 J j-p ^-p m k n G mj (k -p,u- uj')C ni (p,uj') (A8) 

where C ni = G nr CoG* ri . This is the generalization of the self-consistent mode coupling 
of the incompressible turbulent flow, to the compressible turbulent flow and is the mode 
coupling version of the renormalization group arguments of Staroselsky et al. (1990). As 
in all such problems, the mode coupling integral is valid over a larger momentum scale 
and hence, in principle, allows more than the asymptotic analysis of the renormalization 
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group. In this work, the focus is only on the exponent a, which is an asymptotic result. 
In the incompressible limit, it is to be noted that = P^—iuj + vk 2 )~ l , which forces 
Cij = PijC (uj 2 + z> 2 /c 4 ) _1 /c _ ( D+a \ where C is a constant. This reduces equation (A8) to 

~fk ^ — f ^ f k mknP mj {^ - p)Pm(p)P j l(k)P li (k) 

U[ ,U) k 2 J (2ir) D J 2ir [- i(lu - lu>) + v(k- P ) 2 ][lu 2 + t> 2 p*}p D + a 1 ' 

which is very close to the expression obtained by Bhattacharjee (1991). The Kolmogorov 
spectrum for incompressible turbulence requires that a — 0, which characterizes the nature 
of the energy transfer between the mean and the random flow. This transfer characteristic 
should be independent of the speed of sound, and so what holds for c s — > oo, should also 
hold at finite c s . Consequently, the forcing function is characterized by a = in the nearly 
incompressible regime that is being studied here. 

A scaling analysis of equation (A8) can now be carried out. The left hand side scales as 
k 2 ~ y . The right hand side clearly scales as k D+2 - y+2 k y ~ 2 k~ D k 2iy - 2) = k 2y ~ 2 . For the scaling 
properties of the right and the left hand sides to agree, it is necessary to have the condition 
2 — y = 2y — 2, which gives y = 4/3. This leads to the result v ~ A; -4 / 3 , Jx ~ k~^ 3 and 
c s ~ k~ l '\ 
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